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Abstract 

We study the rare decays B+ 7r+£+^", p+i+T and = e,/i) in 

the R-parity violating supersymmetric standard model. From the latest upper limits 
of B{B^ 7r^£~^£~) and B{B^ £~^£~), we have derived new upper bounds on the 
relevant R-parity violating couplings products, which are stronger than the existing ones. 
Using the constrained parameter space, we present the R-parity violating effects on the 
branching ratios and the forward-backward asymmetries of these decays. We find that 
B{B'^ l'^£^) and B{B^ p^£^£^) could be enhanced several orders by the R-parity 
violating sneutrino and squark exchanges, respectively. The R-parity violating effects on 
the dilepton invariant mass spectra of B^ tt'^£'^£~ , p'^£'^£^ and the normalized forward- 
backward asymmetry Afb{B^ — > tt'^£'^£~) and ApsiB^ — > p'^£~^£~) are studied in 
detail. Our results could be used to probe the R-parity violating effects and will correlate 
with searches for the direct R-parity violating signals at the future experiments. 
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1 Introduction 



In the Standard Model (SM), rare B decays B+ ii+l+t, p+i+r and 5)^ ^ i+r {i = 
e,fi) are induced by 6 ^ di~^i~ flavor changing neutral current (FCNC) and expected to 
be highly suppressed. These decays could be an important testing ground of the SM and 
offer a complementary means to search for new physics by probing the indirect effects of new 
interactions. 

The physical aspects of these decays are similar to the b — > si~^i~ decays, which, due to 
its much larger rates, have been studied much more extensively in the SM [H [21 El H] and its 
various extensions [SI El [71 El [9], [10] . Additionally to the b s£~^i~ decays, the b — ^ di~^i~ 
decays would serve an independent test for new flavor changing interactions. Experimental 
progresses aimed to such semileptonic FCNC decays have been made by BABAR [11] and Belle 
[12] with the measurements of the forward-backward asymmetry {Afb) oi B K*i~^i~ and the 
branching ratios of both B —>■ K*i~^i~ and K£~^i~ . Furthermore, the radiative b ^ d'-f penguin 
processes B ^ p'-f and u'-f have been measured by BABAR [13] and Belle [11]. Along with these 
progresses, BABAR, Belle and CDF have begun to probe b di~^i~ decays. Recently, BABAR 
has made a search for B 7r£^£~ and B^ £~^i~ decays [ISj. It's interesting to note that the 
new upper limit of B{B — > iri'^i") < 9.1 x 10"*^ at 90% C.L. [15] has improved the previous 
upper limits [T7j by four orders of magnitude. Meanwhile, the upper limit of B{B'^ f^~^f^~) 
has been improved to be ^3(5)^ l^^l-^") < 1-5 x 10"* at 90% C.L. by CDF [H]. These renewed 
upper limits from BABAR [15] and CDF [16] will provide powerful constraints on new physics. 

In the literature, many useful observables in i? ^ pi~^£~ and B tiI^I~ decays, such as 
branching ratio, CP-asymmetry, forward-backward asymmetry etc., have been investigated in 
the framework of the SM [18], the two Higgs doublet models [19l [20l [21] and supersymmetry 
models (SUSY) [22]. With the new upper limits, it would be very worthy to study these decays 
in other new physics models to derive bounds on the relevant parameters. 

In this paper, we will study these decays in the R-parity violating (RPV) supersymmetric 
standard model [23 l[Ml[25] . The phenomenological constraints on the RPV couplings have been 
studied extensively in B decays [26]. At first, we will update the SM expectations of these decays 
with the up-to-date inputs, such as the new light-cone QCD sum rules results for B 7r(p) 
form factors [27] and the electro- weak parameters [28]. To make realistic estimations, we have 
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taken into account of the uncertainties of the input parameters by varying them randomly 
within la variance. 

The decays 5^ p^£^£~ , vr+^+f" and — ^ £^£~ are all induced at the parton level 

by the same set of the RPV coupling products in the RPV SUSY. From the latest experimental 
data and the theoretical parameters, we will derive the new upper limits on the relevant RPV 
coupling products. Then we will use the constrained parameter space to predict the RPV 
effects on the branching ratios of 5^ — > p^i^i^, B^ 7i~^i~^i~ and B^ decays, 
and the forward-backward asymmetries of B^ p^i'^i" and B^ — > 7T~^i~^i~. Moreover, we 
will compare the SM predictions with the RPV predictions about the dilepton invariant mass 
spectra and the normalized forward-backward asymmetries in B^ p'^i'^i" and B^ n^i^i" 
decays. 

The paper is organized as follows. In Sec. 2, we derive the expressions for B^ —>■ p'^£~^£~, 
B^ 'K~^i~^i~ and B^ £~^i~ processes in the RPV SUSY. In Sec. 3, our numerical analysis 
are presented. We display the constrained parameter spaces which satisfy all the upper limits 
of B{B^ 7r+£"^£~) and B{B^ £'^£~), and then we use the constrained parameter spaces 
to present the RPV effects on Afb{B+ p+^+r), Afb{B^ ^ 7r+£+r), B{B+ vr+^+f"), 
B{B^ p~^£~^£~) and B{B^ t'^t^). We also show the RPV effects on the dilepton invariant 
mass spectra and the normalized forward-backward asymmetries in B^ p+f+f" and B^ 
7i~^£~^i~ decays. Sec. 4 is devoted to our summary. 

2 The theoretical frame for B+ ir+i+i-, p+i+i' 
and B^^ 

2.1 The decay branching ratios in the SM 

2.1.1 The semileptonic decays 5+ n+i+i- and B+ p+E+i- 

In the SM, the effective Hamiltonian for the b di~^£~ decay can be written as [31 [221 ED] 

4^ C 10 2 ^ 

^7/ = -^V^<^Vt, |Ea(p)0,(p) + A:EC,(p)[0,(p) - Ot\f^)] I , (1) 
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where A* = 



The exphcit form of all operators Oi{fi) and the Wilson coefficients Ci{fi) 



calculated in the naive dimensional regularization (NDR) scheme can be found in [3]. The 
effective Hamiltonian leads the QCD corrected matrix element for b d£~^£~ [3], | 



M^^\h dtr) 



V^TT 



(2) 



with Pl,r = (1 =F 75)/2, s = and g = p+ +p_ {p± are the four-momenta of the leptons). We 
take rrid/mh = 0, but keep the leptons mass. The hat denotes normalization in terms of the 
B-meson mass, m^, e.g. s = s/m^, rhg = m^/m^. 

Cg'^'^{fi) contains the one-gluon corrections to Og {uj{s) term) [21] and the corrections of 
four-quark operators Ci_6 and in Eq.(l) [32], which can be written as [29] 



2, 



(3) 



with 



riNDR 
'^9 



UJiS 



g{m,, s) (3Ci + C2 + 3C3 + Q + SC^ + Q 
1 



-(7(1, s)(4C3 + 4^4 + 3C5 + Ce) - -9(0, s)[C, + 3C, 
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g{m,,s)-g{0,s)]{3Ci + C2). 



(4) 
(5) 



In addition to the short distance contributions, b — > d£~^i~ decay also receives long distance 
contributions from both on- and off-shell vector mesons, which can be incorporated in Cg-^'^ by 
the replacement [33l [3] 



g{mc,s) 2^ K 



(6) 



This issue has been extensively discussed in literature [351 EH] . 

With these formulas, one can derive the decay amplitudes for 7T~^i~^i~ and p~^i~^i~ 

decays. Recently the form factors for 5 ^ tt and B ^ p transitions have been improved in 
the framework of light-cone QCD sum rules [27], with one-loop radiative corrections to twist-2 
and twist-3 contributions, and leading order twist-4 corrections. Using the form-factors defined 



in [27] and the matrix in Eq.(l2]), we get the following amphtudes for 

B+ p-^i+e^ decays: 



TT 



M^^\B Mi 
where for Rt vr+f+f", 



2V2 



TT 



and for B+ p+i+r, 



= E{s)e,p^pe*Pp'^sfp-iF{s)el + iG{s){e* ■Pb)P^, + iH{s^^^^ 

with p = pb + Pm{M = 7r+ or p"*" ). 

The auxiliary functions in 7i are defined as [6] 



B\s) 
C\s) 
D'{§) 
A{s) 

B{s) 

C{s) 

m 

E{s) 

m 

G{s) 
H{s) 



1 + 

ci"{s)f.{s) - ^(1 - m^)a;^^Ms), 

Cio/+(S), 

1 + rrin s 



{I + Trip 
1 

1 — 



C'/fis)A,is) + ^(1 - mp)C'/%is) 



(1 - mp)C'/'{s)A2{s) + 2mbC'/' n{s) 



Cl"{s){{l + m,)Ai(s) - (1 - mp)A2{s) - 2mpAo{s)) 
-2mtC'/^T3{s) 
-CioV{s) 



1 + m, 



{l+mp)CioAi{s) 
-^—CioA^is), 



s 



(1 + mp)Ai{s) - (1 - mp)A2{s) - 2mpAo{s) 



The kinematic variables fs, u) are chosen to be 



which are bounded as 



s = f = (p_^+p_y^ 



(24) 
(25) 



{2mif < s < (1 - m 
—u{s) < u < u{s), 



PJ ' 



(26) 
(27) 



with rhi = nii/mB and 



u[s] 



A 



2s-2mlAl + s) 



{2i 



(29) 



Note that the variable u corresponds to ^, the angle between the three-momentum of the 
f"*" lepton and the B meson in the dilepton center-of-mass system (CMS) frame, through the 
relation u = —u{s)cos6 |37j. Keeping the lepton mass, the double differential decay branching 
ratios B"^^ and B^^ for the decays 7T~^i~^i~ and B^ p'^£~^£~ , respectively, are found to 

be 



dsdu 



dsdu 



tb- 



2ii7r5 



■\vuv:,\'{m'+\cr){\-u') 



+ \C'\'Ami{2 + 2<+ - s) + Re{C'D'*)8mi{l - m^+) + iD'l'Amjrs 



(30) 



rB- 



2i%5 



\VuV, 



*|2 
tb\ 



X 



-25^ 

A 



^(s(A + w^) + 4m2A) + ^(5(A + w^) - Amjx) 
f (a - + 8m2+(s + 2mj)) + |Fp(A - ^ 

Re{BE*) + i?e(AF*)" 
[|C|2(A - u^) + |G|2(A - + 4m2(2 + 2mJ+ - i 



^m^+(s — 4m^ 



.^[i?e(SC*)(l-mJ+-s)(A-^2) 



' " 2 

m 



+ Re{FG*) ({I - ml+ - s){\ - u^) + Am]\ 

2 ^2 



Re{FH*) -Re{GH*){l-m'+) +\H\'^sX}. (31) 



V+^J ' I"" I ^2 



2.1.2 The pure leptonic decays — > 

In the SM, the purely leptonic decays B^ — > proceed trough electroweak penguin diagrams 
with Z exchange as well as W box diagrams. The photon penguin contribution is forbidden by 
the conservation of electromagnetic current. The effective Hamiltonian for these decays is [3H] 

where Xt = mf/m'^Y, which leads to the branching ratio 



167rmB^ ^ 

In the SM, as shown by Eq.( l33i) . the decays are suppressed by helicity conservation (the factor 
m1/m\ ) in addition to the small Vtd- However, they are very sensitive to new physics with 
pseudo-scalar interactions. 



2.2 The decay amplitudes in the RPV SUSY 

In the most general superpotential of the minimal supersymmetric standard model, the RPV 
superpotential is given by [23], [Ml ES] 

W^, = ^i:UH^ + \\,^,Ul,EI + K,^UQ,bl + \>:i^^,ptDpl (34) 

where z, j, = 1, 2, 3 are generation indices and c denotes a charge conjugate field. Li{Qi) 
are the lepton(quark) SU{2)l doublet superfields, E'', tJ^ and Z)'^ are the singlet superfields, 
respectively. The bilinear RPV superpotential terms fiiLiHu can be rotated away by suitable 
redefining the lepton and Higgs superfields [25]. However, the rotation will generate a soft 
SUSY breaking bilinear term which would affect our calculation through penguin level. The 
processes discussed in this paper could be induced by the tree-level RPV couplings, so that we 
would neglect sub-leading RPV penguin contributions in this study. The A and A' couplings 
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in Eg. (1341) break the lepton number, while the A" couphngs break the baryon number. There 
are 27 A'^^^ couphngs, 9 Xijk and 9 Af^^ couphngs. X[ij]k are antisymmetric with respect to their 
first two indices, and X'-y/^^ are antisymmetric with j and k. All the processes considered in 
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Figure 1: The RPV contributions to 7r~^ {p~^)i'^i~ due to the sneutrino and squark 

exchange. 





(a) 

Figure 2: The RPV contributions to -Bj] 



(b) 



due to the sneutrino and squark exchange. 



this paper involve terms in the RPV Hamiltonian with two leptons and two quarks as external 
fields. From Eq. flM|) . the Hamiltonian of 6 ^ d£^£~ process due to the squarks and sneutrinos 
exchange are 



^eff 



1 ^ ^2ih^^d,YPRdn)iin,PLi,] 



E 



2 {d^PBdn)iikPLij) + ^^^^^{dnPLdmmPRik) 



m 



(35) 



The RPV Feynman diagrams of ir^i^i^ , B^ p~^i~^i~ and B'^ — > are displayed 

in Fig{T]and Figl21 respectively. 

From Eq.( !35l) . we can obtain the RPV decay amplitude for B^ tt"' 



M^^ {Bl 



TT 



'k{PB+ Pn+ 



)(1 - 75)^,) + T,(4(l - 75)^,) + n(4(l + 75)4),(36) 
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where 



T- 

J- 7 



5m: 



4m? 



m. 



2 



nib - rrid 



\* \/ 



m 



B 



2 

mf.+ 



mt - md 



with mfo and m^ are the quark's running masses at the scale m?,. 
For 5+ ^ p+£+£^, the RPV amphtude is 



(5+ ^ p+£+r) = r3^(47^(i - 75)^i) + *p(4(i - 75)^,) + $',(4(1 + 75)^,), 



where 



E 



E 



2 m^ + m^ 



-A2m| 



A^m^ 



2 mfe + m^ 

and the auxihary functions above are found to be 



m 
m 

K{s) 

m 



E 



5m,- 



niB + m. 



-m. 



E 



5m, 



(ms + mp+)A^ 



B,, 



E 
E 



5m: 



5m: 



2m^ 



-m 



B 



For the pure leptonic decays 



the RPV amphtude is 



M"^^ {Bd ^ 44) = fi^(4 - 75)4) + ^^(4(1 - 75)4) + ^'u{W + 75)4). 



where 



^/Sd^Bd E 



% = ^fB.f^B, E ^ 
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Generally, the RPV couplings can be complex and their phase may induce new contributions, 
so we write them as 



and G [— vr, vr] is the RPV weak phase. 



\* \ I \ \* I p—ivifp 



(50) 



2.3 The branching ratios with RPV contributions 

Using the formulae presented above, we get the following results for the double differential 
decay branching ratios for the decays M~^i~^i~ (M"*" = n'^ or p^), 

d^B^'^ d^B'£ d^B¥^ d^B',^'^ 



(51) 



dsdu dsdu dsdu dsdu dsdu 

Since it is always assumed in the literature for numerical display that only one sfermion con- 
tributes at one time, we have neglected the interferences between different RPV coupling prod- 
ucts, but kept their interferences with the SM amplitude, as shown in the following equations. 

For the 5+ 7r+£+£- decay, 

d'^Bf 



ml 



dsdu 



Re{WA'Tl){X~u' 



277r3 { 

+ReiWC'Tl) [ - (A - n^) - Am^,{2 + 2ml+ - S) 
+Re{WD'T*-) [ - 4m2(l - ml+) 

+ |TspmB [a - + 2mj{2 + 2ml+ - s] 



(52) 



rrii 



dsdu 



•2 KS'n^ 



277r3 



Re{WA'Tl){2meu) + Re{WC'Tl){l - ml+){-2mi) 



+Re{WD'Tl){-2mis) + |T^P(s - 2m 



I) 



(53) 



d'B\ 



dsdu 



Tb 



^B 

277r3 



Re{WA'T;){2miu) + Re{WC'T;){l - ml+){2me) 



+Re{WD'T';){2mes) + |T'-|^(s - 2m 



(54) 



and W = -^VlVurriB. 
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For the B+ p+i+i' decay, 



dsdu 



d^Bf 
dsdu 



d'^B't 
dsdu 



Re{WAr) s(A + n^) + 4mj\ 



-ReiWEF) [s{\ + n^) - Am^^X] + |/p [s(A + n^)' 

+Asu\Re{WAJ*) + ReiWBF) - Re{WEJ*) - i?e(VrFr) + 2i?e(JJ* 



m: 



+ 1 J|2(A - + 8mJ+(s - m|)) 
-Re{WBK*){X - u^){l - ml+ 



+Re{WFK*)(^{\ - u^) (1 - ml+ - s) + Am^^X 

-2Re{JK*)({X - u^){l - ml+ - s) + 2m|A) 
A 



m: 



Re{WCK*){X - r) - Re{WGK*){X - u' + 4m|(2 + 2m 



+i^r(A 



2m^(2 + 2m 



4m» 



m: 



Re{WHL*)s + \L\'s/2 + Re(WFL*) - Re{JF 



-Re{WGF){l - ml+) + Re{KL*){l - ml 



m. 



ml 



+Im{WC^l)X^u + Im{WF^)X 



Im{WB^) ■ ( A-2m(1 -ml+-s 



-Im{WG<^l)X^{l-m^ 



Pfs-2m2) 



Tb 



277r3' 



m; 



m: 



Im{WB<^>';) ■ (X^h{l - ml+ - s) 



+/m(iyC$'-*)A2n - Im{WF<l>';)X'' 



+Im{WG<i>';)X^{l-m 



2mi 



The normahzed forward-backward asymmetries is defined by 
Afb{B Mi+r) = J ds 



I-i '"'^^f^cr/'"^ stgnicose)dcose 
1+1 fm^^lM^rfco.^ 



dsdcosO 
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Since there is no term containing u with an odd power as shown by Eq.( l30i) . the Afb is zero 
for TT~^i~^£" decays in the SM. The RPV effect via squark exchange on Afb{Bu 

7r~^ i'^i") also vanishes for the same reason as shown by Eq. fl52p . while the sneutrino exchange 
contributions to Afb, as shown by Eqs. (l53l 15^ . are proportional to m^, which are tiny for 
i = e, fi. 

The total decay branching ratios of the pure leptonic decays are given by 

„, „ . f 1 



B''\B, 
1 

1 



\KsM 

1 



I a 



mi. 



V2m^ 



\K 



with 



K. 



sm\ 



-Re{KsM%*)— + \^,\ 
mi 



,/ 1 



\2m} 



m 



a. 



SM 



-VuV:,Y{xt){ifB, 



(59) 



(60) 



y/2 2TTsm29w 

From Eq. fl59l) . we can see that B{B^ £^£^) could be enhanced very much by the s-channel 
RPV sneutrino exchange, but not so much by the t-channel squark exchange. 



3 Numerical results and analyses 

With the formulae presented in previous section, we are ready to perform our numerical analysis. 
Firstly, we will specify the input parameters. Then we will use the latest up-limits of B{B^ 
7r~^ and B{Bj^ i'^i^) to get the constraints on the relevant RPV couplings. Finally, we 
will show the distributions of the branching ratios and the forward-backward asymmetries of 
these decays in the constrained RPV parameter space. 

In the Wolfenstein parametrization of the CKM matrix elements, one has 

K-^-j^,^oin ,61) 

The parameters X,p and f/ have been updated by CKMfitter Group [39], which read 

A = 0.22568^™, -p = 0.213^°:°^t v = 0.3481°:°^?. (62) 
In our numerical analysis, we have used IV^^^V^^I — \ Vtd\ — (8-27l^g'49 
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For the semileptonic decays ix^l^l" and — > p~^i~^i^ , we will use the form factors 

of light-cone QCD sum rules (LCSRs) results ^7\, which are renewed recently with radiative 
corrections to the leading twist wave functions and SU(3) breaking effects. For the depen- 
dence of the form factors, they are parameterized in terms of simple formulae with two or three 
parameters. The form factors V, Aq, Ti, /_|_ and fx are parameterized by [27] 

For the form factors A2 and T3, it is more appropriate to expand them to the second order 
around the pole 

^^'^ ^ l-s}m%,^ il-s/m%r ^^^^ 
The fit formula for Ai, T2 and /o is 

^ T^,- 

The form factor T3 can be obtained with T^ls) = ^—p^[T3{s) — T2{s)]. The values of all the 
corresponding parameters for these form factors are listed in Table I, the uncertainties induced 
by F(0) are also considered in the following numerical data analyses. 

The values of other input parameters used in the numerical calculation and the experimental 
upper-limits are given, respectively, in Tables II and III. In the numerical calculations through 
the paper, we will take into account of the la uncertainties of the parameters with relative 
large error bars. 

With these inputs, we get the SM prediction B{B+ 7r+£+£-) = (2.03 ± 0.23) x 10-^ 
which is smaller than the result (3.27 x 10~^) in Ref. [I9]. The difference is mainly due to the 
updated values for formfactors and CKM parameters. 

As in the literature, we assume that only one sfermion contributes at one time with a mass of 
100 GeV. For other values of the sfermion masses, the bounds on the couplings in this paper can 
be easily obtained by scaling them by factor = (lo^fcy)^' present the branching ratios 
of the SM expectation in Table III, where the relevant experimental upper limits [inillSlSI] are 
listed for comparison. We can see that the experimental upper limits of B{B^ TT~^i~^i~) are 
just one order above the SM expectations, which will constrain RPV SUSY parameter space 
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Table I: Fit for form factors involving the B 7r(p) transitions valid for general q"^ |27j . 



F{s 


) 


m 


Atot 


ri 


m| 


r2 




fit Eq. 




■TV 


0.258 


0.031 


0.744 


5.322 


-0.486 


40.73 


m 




■IT 


0.253 


0.028 


1.387 


5.322 


-1.134 


32.22 


m 




■TV 


0.258 


0.031 






0.258 


33.81 




yB^ 


*P 


0.323 


0.030 


1.045 


5.32^ 


-0.721 


38.34 




At 


>P 


0.303 


0.029 


1.527 


5.28^ 


-1.220 


33.36 


m 




*P 


0.242 


0.023 






0.240 


37.51 


m 




♦P 


0.221 


0.023 


0.009 




0.212 


40.82 




rjnB- 


♦P 


0.267 


0.023 


0.897 


5.322 


-0.629 


38.04 


m 


Tt 


♦P 


0.267 


0.023 






0.267 


38.59 


m 


ft 


♦P 


0.267 


0.024 


0.022 




0.246 


40.88 


dMD 



Table II: Default values of the input parameters and the itlcj error bars of the sensitive 
parameters used in our numerical calculations in addition to the ones discussed in text. 

niB^ = 5.279 GeV, tub^ = 5.279 GeV, mw = 80.425 GeV, 

= 0.140 GeV, m^o = 0.135 GeV, nip = 0.776 GeV, 
mbirrib) = (4.20 ± 0.07) GeV, 

mu{2GeV) = 0.0015 ~ 0.003 GeV, md{2GeV) = 0.003 ~ 0.007 GeV, 
me = 0.511 X 10^3 Q^y^ ^ 0.106 GeV, rut^poie = 174.2 ± 3.3 GeV [2H] 
tb^ = 1.530 ps, tb^ = 1.638 ps. [2H] 
sin^Ow = 0.22306, = 1/137. [2H] 
fB, = 0.2161°:°?^ GeV m 



significantly. Although the present limit of B{B^ /^^/^~) are three orders of magnitudes 
higher than the SM predictions, the upper-limit will also yield strong bounds since the decay 
fi^fi" is very sensitive to new pseudo-scalar interactions beyond the SM. 
Now we turn to the RPV effects. The decays to fi+ p+i+i', fi+ vr+f+f" and 
^ i+i- involve the same set of the six RPV coupling products. Combining the upper- 
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Table III: The SM predictions and experimental data of branching ratios for the 
upper limits for £+£- and B+ tt+£+£- [T^fTniliT]. 





SM prediction value 


Experimental data 






) 


(2.03 ±0.23) X 10^^ 


< 2.8 X 10-^ (90% C.L.) 


m 


mt - 


-> vr 


^e+e") 


(2.03 ±0.23) X lO^s 


< 1.8 X 10-^ (90% C.L.) 


m 


mt - 


P 




(4.33 ± 1.14) X lO^s 






mt - 


P 




(5.26 ± 1.37) X 10^^ 






B{B,- 






(1.33 ±0.12) X 10-1° 


< 1.5 X 10-8 (90% C.L.) 


m 


B{B,- 






(3.11 ±0.27) X 10-15 


< 6.1 X 10-8 (90% C.L.) 


m 



5 
> 

Q. - 
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Figure 3: The allowed parameter spaces for the relevant RPV couplings constrained by the 
measurements listed in Table IV, and the RPV weak phase is given in degree. 
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limits of B{B+ ^ iT+i+e~) p] and i+E^) [lElEH] listed in Table III, we get constraints 

on the six RPV coupling products. Subject to the constraints, the random variations of the 
parameters as discussed previously lead to the scatter plots displayed in FigOl From the figure, 
we can see that the constraints RPV on Ag^aAg*! and A'^^gA']*^ are sensitive to their RPV weak 
phases, which implies significant interferences between the squark contributions and the SM 
ones. As depicted by Fig.l and 2, the squark contribution amplitudes are of b{V + A)di{V — A)£ 
current structure after Fierz transformations, which are similar to the SM amplitudes. For the 
constraints on A'A* of sneutrinos are insensitive to their phase . We summarize all the 
constraints in Table IV, where the recent updated bounds [12] are listed for comparison. It's 
shown that our bounds are more restricted than the previous ones. We note that the strong 
constraints on yu^Xin and A2j3A2*i are mainly due to the upper-limits of B{B^ n^i^i^), 
which are just one order higher that the SM expectations, and the constraints on other four 
RPV coupling products from the sneutrino exchange are mainly due to the upper-limits of 

Table IV: Bounds for the relevant RPV couplings products by ir^l'^i" 
and B'^ — > i~^£~ decays for 100 GeV sfermions, and previous bounds are listed for 
comparison. 



Couplings 


Bounds [Processes] 


Previous bounds 


Processes] |l42j 


l^liS^lill 


< 1.1 X 10-^ 


[Bd- 
Wt- 


■*e+e-] 
->7r+e+e-] 


< 2.6 X 10-2 [Bd 


^ e+e-] 


AjiiA-31 


< 7.4 X 10-6 


[Bd- 
[Bt- 


■+e+e-] 
->7r+e+e-] 


< 4.1 X 10"^ [Bd 


^ e+e-] 


1 \* V 1 

l^ill^ilSl 


< 7.4 X 10-6 


[Bd- 

[Bi- 




< 4.1 X 10-^ [Bd 


^ e+e-] 


1 V \'* 1 


< 1.3 X 10-^ 


[Bd- 

[Bt' 




< 5.4 X 10-^ [Bd 


^ /i+/i-] 


Ai22A-3i 


< 3.7 X 10-6 


[Bd- 
[Bt' 




< 6.2 X 10-6 [Bd 


/i+/i-] 


I^j22^il3l 


< 3.7 X 10-6 


[Bd- 
\Bt' 


-*TT + fj,+ H^] 


< 6.2 X 10-6 [Bd 


^ /i+/i-] 



Using the constrained parameter spaces shown in Table IV and Fig|3l one can predict the 
RPV effects on the other quantities which have not been well measured yet in these processes. 
We perform a scan over the input parameters and the new constrained RPV coupling spaces 
to get the allowed ranges for B{B^ P+^+^~), B{B^ 7r+£+£-), B{B^ ^+r), ApBiB^ ^ 
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7i~^i~^i ) and Afb{B:^ — ^ p^i'^i ) with the different RPV couphngs products. The numerical 
results are summarized in Table V. 

Tabfe V: The results for B{B^ i+£-), B{B+ p+e+£- ,Tr+e+£-) and Afb{B+ 

in the SM and the RPV SUSY. The RPV SUSY predictions are obtained by the allowed regions 

of the different RPV couplings products. In the RPV coulpings, g = I and 2 for ^ = e and p, 

respectively. 





OlVi VdlLlL- 




\ X'* 


X* \' 

'^igg-^ilS 


B(B" e+e-) 


[2.2,3.7] X 10-15 

L ' J 


[0.002,6.7] X 10-" 


[0.008,6.1] X 10-* 


[0.008,6.1] X 10-* 


B{B+ TT+e+e-) 


[1.4,2.6] X 10-8 


[0.57,18.0] X 10-* 


[1.6,15.5] X 10-* 


[1.7,13.8] X 10-* 


B{B+ ^ p+e+e-) 


[2.8,9.7] X 10-* 


[3.8,586.8] X 10-* 


[2.9,10.2] X 10-* 


[2.8,10.2] X 10-* 


Afb{B+ ^ TT+e+e-) 








[-7.2,6.7] X 10--5 


[-7.2,7.2] X 10-5 


ApBiBi ^ p+e+e-) 


[0.14,0.23] 


[-0.13,0.82] 


[0.14,0.23] 


[0.15,0.23] 




[0.96,1.56] X 10-1° 


[0.01,38.9] X 10-1° 


[0.01,1.5] X 10-* 


[0.01,1.5] X 10-* 


B{B+ ^ n+p+p-) 


[1.4,2.6] X 10-*^ 


[0.22,27.92] X IQ-* 


[1.7,5.3] X 10-* 


[1.6,5.4] X 10-* 


B{B+ ^ p+p+p-) 


[2.8,8.4] X 10-* 


[3.5,780.4] X 10-* 


[2.8,9.1] X 10-* 


[2.8,9.1] X 10-* 


AFBiB+ TT+p+p-) 








[-1.1,1.1] X 10-2 


[-1.1,1.1] X 10-2 


AFsiBi ^ p+p+p-) 


[0.14,0.23] 


[-0.13,0.46] 


[0.14,0.23] 


[0.14,0.24] 



For Table V, we have the following remarks. 

• As shown by Figl2]J^b), the contributions of ^1^3X^1 and A2j3A2*i to B{B^ e+e~) and 
B{B^ —>■ fi+p,~), respectively, arise from t-channel squark exchange. After Fierz transfor- 
mation, the effective Hamiltonian due to the t-channel squrk exchange is proportional to 
b'^^P^di'-y^PLi, which contributions to B{B^ — > i+i~) are suppressed by m1/m\ due to 
helicity suppression. Therefore, B{B^ — ^ e~^e~ , fi~^ p.~) will not be enhanced so much by 
the t-channel squark exchanging RPV contributions. In other words, the bounds on these 
two RPV coupling products are due to the upper-limits of B{B+ — > n+i+i'). With the 
bounds, the RPV squark exchanges could enhance B{B'^ e+e~) and B{B^ p+fi"), to 
the utmost, to 6.7 x IQ-^^ and 3.9 x IQ-^, respectively. However, the effective Hamiltonian 
of s-channel sneutrino exchange would be b{l ± '~f^)d£{l =F 75)^, whose contributions are 
not suppressed by ml/m^. Therefore, both B[B^ — > e+e-) and B{B^ 1^^ could 
be enhanced to order of IQ-^ to saturate their experimental upper limits by XniX^i and 
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Figure 4: The effects of RPV coupling X'li-^Xin due to the squark exchange in n'^e'^e' , 

— s> p^e^e~ and B^ e~^e~ decays. The primed observables are given with ip^nS) VMD 
contributions. 

Ai22A-3i, respectively. 

• At present there is no experimental data for B^ p^i'^i^ decays. It is interesting to 
note that B{B^ p~^i~^i~) could be enhanced to order of 10~^ by the t-channel squark 
exchange within the bounds due to the upper-limits of B{B^ 7i^i^i~), however, it 
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Figure 5: The effects of RPV coupling \'2iz'^'2ii due to tlie squark exchange in n'^fi^fi , 

— * p'^fi'^fi" and B^ — > decays. 

is not sensitive to the s-channel sneutrino exchanges within the constrained parameter 
space. 

• Afb{B:^ -n^e^e~) and Afb{,B^ n^fi^fi") are zero in the SM. The RPV contri- 
butions to the asymmetries due to squark exchanging are also zero, while the sneutino 
exchanging RPV contributions are also too small to be accessible at the LHC 
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• It is found that the squark exchanging RPV contributions have significant impacts on 
Afb{B^ p"^e^e~) and Afb{B^ P^/^^/^")- While, the RPV contributions due to 
the sneutrino exchange affect Afb{B:^ ^ p^e'^e") and Afb{B:^ p^fi~^fi^) shghtly. 

For each RPV couphng product, we present the distributions and correlations of branching 
ratios and forward-backward asymmetries within the constrained parameter space in FiglHl 
Numerical results are shown in FigsHHHl where correlations between the physical observable B, 
Afb and the parameter spaces of the different RPV coupling products by the three-dimensional 
scatter plots. The dilepton invariant mass distribution and the normalized forward-backward 
asymmetry are given with vector meson dominance (VMD) contribution excluded in terms of 
dB/ds and dAFB/ds, and included in dB'/ds and dA'FB/ds, respectively. From FigsJlMl one 
can find the correlations of these physical observable with the RPV coupling products. 

Now we turn to discuss plots of FigH] in detail. The three-dimensional scatter plot shows 
Afb and B correlated with |A'^j3A']*^| and its phase (p^^. We also give projections on three 
vertical planes, where the lA'^^gA'^* J-^^^ plane displays the constrained regions of A'^^gA'^*]^ as the 
second plot of FigI31 From plots of Figjl](a)-(c), we can find that B{B^ e"''e~), B{B^ 
Ti^e^e~) and B{B^ — > p~^e~^e~) are increasing with lA'j^jgA']*^ on the B-lX'^-^Xi^^l planes, and 
B{B^ 7r"^e~^e~) can be enhanced to its present experimental upper limit. Fig|4](d) shows that 
Afb{B:^ — i> p+e+e~) is decreasing with |A'^j3A'^*]^| on the Afb{B:^ P^^~^^~)'\Ki3'^iii\ plane. 
It is shown that B{B'^ e"'"e~), B{B^^ p^e^e^) and Afb{B:^ p~^e~^e~) are sensitive to (p^^ 
with the first two increasing with \(f)^^ \ while the last one decreasing, whereas B{B^ Ti^e^e~) 
is insensitive to (f)^^. From FigsJll^e)-(f), we can find that the the A'^^^gA'^*^ contributions are 
significantly larger than the SM contributions at the all s regions with theoretical uncertainties 
included. While for dB{B^ p'^e^e^)/ds, the RPV effect is significant at the low s region. 
Especially from Fig|4](i)-(j), we find the RPV contributions to dA^pBiB^ p^e^e~)/ds are 
quite different from the SM expectations. All the above observations can be extended to FigEl 
but for squark UiL contributions to B^ P^/^~) 7r^/i"^p^ and B^ p^p~^p~ , instead. 

The RPV sneutrino exchange contributions to i?^ p^£^£^, B^ 7T~^i~^i~ and B^ 
i~^£~ are very similar to each other. We would take A*22A^]^3 contributions as an example, which 
is shown by FigJHl Figsin](a)-(b) show that both B{B^ /^^P~) and B{B^ 7i~^fi~^fi~) are in- 
creasing with |A*22A^]^3| on the 'B-|A*22A^i3| planes, but insensitive to , and B{B^ P^p ) can 
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Figure 6: The effects of RPV coupling A*22A'j]^3 due to the sneutrino exchange in — > 7r+/i+/i , 
p^fi~^fj,^ and — > decays. 

be enhanced to ~ 1.5 x 10"*^ which saturates its present upper hmit. As shown by FigslHI^c) and 
(d), both Afb{B^ — i> p^fi^fi^) and B{B^ p^fi^fi~) are insensitive to the RPV sneutrino 
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contribution. The scatters of points in FigsEJ^c) and (d) are due to theorectical uncertain- 
ties. The similar situation is also found for the distributions of dApsiB^ — > fi^ fi~)/ds 
and dB{B^ p^i^i~^fi~)/ds. From Figsin](g)-(j), one can find that the present theoretical un- 
certainties in our numerical calculation are still very large, and the RPV sneutrino exchange 
contributions are indistinguishable from the uncertainties. Figsj6](e)-(f) are the distributions of 
dB^'\B^ 7r+yU+/i~)/(is, which show the RPV contributions distinguishable only in the high 
s region. 

As known, to probe new physics effects, it would be very useful to measure correlative 
observables, for example, ^ p+^+r), B{B+ tx+I+I^) and B{Bd since 

correlations among these observables could provide very strict bound on new physics models. 

It should be noted that the present upper limit of B{B^ n^i'^i") < 9.1 x 10"*^ made 
by BABAR [15] is based on just 208.9/6"^ data, which is already just one order higher than 
the SM expectation. Up to this Fall, Belle and BABAR have recorded as much as 709fb~^ 
and 500/6"^ data, respectively. Searches with these data will surely improve the upper limits 
for these decays to large extent, which will give much stronger constraints on the RPV SUSY 
parameter spaces. On the other hand, the decays with enhancement of the RPV SUSY might be 
accessible at BABAR and Belle, especially, the B^ — > p~^£~^£~ decay. For such a measurement, 
a good Tr+Tr*^ sub-vertex reconstruction is necessary to distinguish it from the much more 

copious decay B^ K*~^i~^i~ —>■ {K~^7T^)i~^i~ , and may be a challenging task at BABAR and 
Belle. In the near future at the LHCb, these decays might be measured with much larger data 
samples. 

4 Summary 

We have studied 5^ 7r~^ B^ p~^i~^i~ and B^ — > £~^i~ decays extensively in the 
RPV SUSY. From the latest experimental upper limits of B{B^ 'Ji~^£~^£~) and B{B'^ 

we have obtained the constrained parameter spaces of the RPV coupling products, 
and found these constraints are robust and stronger than the existing ones, which may be 
useful for further studies of the RPV SUSY phenomenology. Using the constrained parameter 
space, we have given the RPV impacts on the branching ratios B{Bj^ — > B{B^ 
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TT^l^l ) and B{B^ p^i^i ), and the forward-backward asymmetries ApsiB^ tt ) 
and Afb{B+ ^ p-^£-^£-). 

It is shown that B{B^ — > is very sensitive to the RPV s-channel sneutrino exchange, 

but not to t-channel squark exchange, which could be enhanced to ~ 10~^ by the former. On the 
contrary, p'^i~^i~ is very sensitive to t-channel squark exchange, and B{B^ p'^i'^i~) 

could be enhanced to ~ 10~^, which could be tested at BABAR and Belle. 

Furthermore, we have shown the RPV effects on dilepton invariant mass spectra and the 
normalized forward-backward asymmetries in B^ n'^i^i' and B^ p~^i~^i^ decays. In 
the constrained parameter space, it is found that the RPV couplings A'^^gA'^*]^ and X'2i3^2ii of 
squark exchanges could result in the distributions of dB{B^ p^i^i~)/ds and dApBiB^ 
p^i~^i~)/ds significantly different from the SM expectations. The other four RPV couplings 
due to the sneutrino exchanges give distinguishable contributions to dB{B^ — > 7r+£+£~)/(is 
at high s region. But for the B^ p^i^i~ decays, the sneutrino exchange contributions to 
dB/ds and dAps/ds are indistinguishable from the theoretical uncertainties. 

With the operation of B factories, many measurements of rare B decays have put strong 
constraints on new physics scenarios beyond the SM. With the near future experiments at 
the LHCb, one can access many more rare decays with high statistics, which will give more 
stringent bounds on the products of the RPV couplings. From measurement of the correlated 
decays, one may eventually establish indirect evidence for specified new physics model or rule 
it out. 

The results in this paper could be useful for probing the RPV SUSY effects and will correlate 
strongly with searches for the direct RPV signals at future experiments. Negative experimental 
evidence for the deviations from the SM expectations in these decays would result in strong 
constraints on these RPV couplings. 
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